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Abstract 

We discuss a general procedure for arriving at the Hamilton-Jacobi 
equation of second-class constrained systems, and illustrate it in terms 
of a number of examples by explicitely obtaining the respective Hamil- 
ton principal function, and verifying that it leads to the correct solu- 
tion to the Euler-Lagrange equations. 
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1 Introduction 



There has recently been much interest in a Hamilton- Jacobi (HJ) formula- 
tion for constrained systems [l]-[3]. Although there exists consensus in the 
formulation for purely first class systems, the case of second-class systems 
has to be dealt with separately. 

For an unconstrained system described by a Lagrangian L(qi,qi), i = 
1, • • • , n, the HJ-equations read, 

H' =p + H (q ,q i ,p i ) = Q (1) 

where H is the canonical Hamiltonian, and where, for later convenience, po 
and pi are shorthand for the partial derivatives 

dS dS 

Po = tt- , Pi = tt- , Qo = t , (2) 
dq dqi 

with S the Hamilton principal function S = S(t, qi) to be obtained as a 
solution of eq. (1). 

Consider now a constrained system in which not all canonical variables 
are independent. In that case, the Lagrangian L(?j, qi) is singular, so that the 
determinant of the Hessian matrix Ha = a 9 ~. is zero and the accelerations 
of some variables qi are not uniquely determined by the positions and the 
velocities at a given time. 

Let n — mi be the rank of the Hessian. In that case only n — mi velocities 
q a (a = mi + 1, - ri) can be solved for as a function of the coordinates 
the momenta p a canonically conjugate to q a , and the remaining velocities 
q a , to yield q a = q a (qi,Pb,qp)- The remaining momenta p a (a = l,...,mi) 
can be shown to be functions of and p a only, and represent the primary 
constraints (f> a (q,p) = Pa + H a (qi,p a ) = 0, in the Dirac terminology [5]. We 
write them in a form appropriate for our purposes: 

H ' a = Pa + H a (qo,qi,p a ) =0, (a = l,---m 1 ) (3) 

Adjoining these constraints to (1), we are led to consider the coupled set of 
differential equations 

dS \ 

Qo, Qa, q a ,^~ = , 

Qo,Qa,Qa, Tj^-J = 0, a = l,...,mi, (4) 
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where Hq is the canonical Hamiltonian evaluated on the subspace of the 
primary constraints, which only depends on the variables g«, p a (and possi- 
bly Qo — t), as is well known. Continuing the above line of reasoning, the 
system of differential equations (4) should further be suplemented with the 
corresponding set of equations associated with possible secondary constraints 
<^(9)P)cr = 0. Let there be m 2 such constraints. Adjoining these to (4), we 
have in our generalized notation, 

dS 

H'J^qi, —) = , i = 0, l,---,n, a = 0, 1, • • • , m = m x + m 2 . (5) 

Assuming that the secondary constraints do not represent constraints among 
the coordinates themselves, 3 this coupled set of differential equations only 
admits a solution provided the H' a are in weak involution in the sense 4 

{^m P ^-CjH' x {q,-), (6) 

where the (generalized) Poisson bracket is defined by 

™ (dAdB 8BdA\ 

Hence, for second class systems, our Ansatz (5) for the HJ equations must 
be incorrect. 

One way of circumventing the integrability problem would be to embed 
the model in question into a larger phase space, such as to become first class 
[8]. Then the integrability conditions (6) are statisfied and our Ansatz is 
self-consistent [6]. This is, however, not the point of view we wish to take. 
As we shall rather illustrate in terms of a number of examples, a direct 
way of proceeding in the case of systems with second class constraints is to 
perform a canonical change of variables in which the second-class constraints 
become part of the variables. If the Poisson brackets of constraints are given 
by constant matrices, one can always perform a transformation to variables, 
where the constraints themselves can be grouped into canonically conjugate 



3 As we shall see, this condition is not satisfied for the "non-linear sigma model" . 
4 On the space of solutions to (4) this then implies strong involution [6]. 
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pairs. Otherwise this may not be possible [7], and the embedding procedure 
would provide a way out of the dilemma. 

Unlike the discussions found in the literature [1, 9, 12] it is the final 
objective of this paper to actually construct the Hamilton-principal function 
S for the examples we consider, and to obtain from there the solution of the 
corresponding Euler-Lagrange equations following standard methods. 

Our strategy to be illustrated by a number of examples will be the fol- 
lowing: We begin by seeking all constraints, following the Dirac algorithm. 
We then seek a canonical transformation from a set of phase space vari- 
ables qi,Pi to a new set (Q,P) = {<ll,Xa,Pl i ,'Pa), in which the complete set 
of primary and secondary constraints are grouped into canonically conju- 
gate pairs (xai'Pa) satisfying {Xa,Vp} = 5 a p- Since the number of second 
class constraints is even, this is always possible for linear constraints. 5 
The remaining set of variables which in turn are canonically con- 

jugate to each other, {q*,pl} = S a b, are chosen to commute with the set 
{XaiVa}- Denoting by F2(q;p*,V;t) the generating function for this canon- 
ical transformation, the new Hamiltonian will be given by H(q*,p*; X , V; t) = 
H(q(q*, p\ X , V),p(q*, p\ X: V);t)+^ : where q(q*, p\ X , V) and p{q\ p\ X , V) 
are obtained by solving the coupled set of equations pi = Xa = §pr and 
q* a = The extended action in terms of the new variables now reads, 6 

S E = J dt{p:q: + V aXa -H(q*,P*,X,V;t)-\ a Xa-VaV a } (8) 

Since the constraints commute with q* and p*, the equations of motion for 
these variables are determined by H with the constraints set equal to zero, 
strongly: 

& = {q* a ,H(q*,P*,0,0;t)}, pi = {p* a ,H(q*,p*,0,0;t)} 

while the equations for the remaining coordinates, the constraints, just repro- 
duce the persistency equations, and represent a set of differential equations 

5 As shown in ref. [7] this is not always possible, if the constraints are nonlinear. 

6 To formulate the general problem it is convenient to work with the extended action, 
where all of the constraints arc implemented via Lagrange multipliers. Since all the con- 
straints are supposed to be generated from the total Hamiltonian Ht, involving only the 
primary constraints, via the Dirac self-consistent algorithm, the Lagrange multipliers as- 
sociated with the secondary constraints will eventually turn out to be zero in a variational 
calculation. 
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in the new coordinates. With respect to the coordinates {q*,p*) the problem 
has thus been reduced to that of an unconstrained system, and the usual 
Hamilton- Jacobi treatment applies. Note that the canonical change in phase 
space variables has now removed the inconsisteny of the originally envisaged 
HJ formulation. The reason is that, unlike the Hamiltonian, the Hamilton 
Principal function in the new coordinates (Q, P) cannot be obtained by a 
simple substitution (q,p) — > (q(Q, P),p(Q, P)). 

In different terms: from the total action (8) we obtain for the HJ equation 

~ , , dS OS . OS n 

together with the constraint equations, 

X« = 0, V a = — = 0. (10) 
In general our problem thus reduces to solving the differential equations 

Atf,g, t) + f-o (ii) 

with 

H(q*,p*) = H(q*,p*, X = V = 0). (12) 

In the following section we illustrate these ideas in terms of a number of 
simple examples. We shall be rather detailed in the first example, in order 
to illustrate various aspects of the problem. 



2 Second class systems and HJ formulation 

In this section we shall consider several simple quantum mechanical models 
of purely second class systems and show how to properly deal with them in 
order to construct the corresponding Hamilton Principal function S. We also 
verify that the formalism yields the solutions to the corresponding equations 
of motion. 
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Example 1 



Consider the Lagrangean 

L = ^x 2 + xy - ^(x -yf (13) 

The corresponding Euler-Lagrange equations can be written in the form, 
x = , x + x — y. We have one primary constraint <fi := p y = 0, and the 
canonical Hamiltonian on the primary surface reads, 

H G = l -{ Px -yf+ l -{x-yf. (14) 
The usual Dirac algorithm yields a secondary constraint f = 0, 

{0,#o} = -2<P, V = V-^(Px + x). (15) 

There are no further constraints. Since {tp, 0} = 1 these constraints are 
second class and canonically conjugate to each other. 

If we were to proceed naivly and make the replacements p y — ^ in the 
primary constraint p y = 0, then this would imply that S is independent of 
y, which would be in direct conflict with the equation obtained by making 
the substitution p x = |^ in the constraint = 0. This reflects the non- 
commutative structure of our second class system, as well as the violation of 
the integrability condition (6). 

In order to circumvent this difficulty, we make a canonical transformation 
to a new set of variables (q*,p*) and (x^X i n which one of the canonically 
conjugate pairs are chosen to be the constraints themselves: 

X = <P, V = 0, 

Q* = x--p y , p* =p x + -p y . (16) 

We evidently have 

{ X ,V} = 1, {q*,p*} = l. (17) 

Futhermore, {q*,p*) and (%, V) "commute". Our change of variables is thus 
canonical. Indeed one readily checks that 

Pxi + p y y = p*q +Vx + \j t (P* v ~ \ v<2 ) ■ 
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Furthermore, the corresponding generating functional F2(x,y;p*,V) for the 
transformation (16) is readily constructed: 

F 2 = (p* - l -V)x + Vy- ^Vp* + ^P 2 . 

The inverse transformation to (16) reads 

x = q* + ^P, p x =P*-^P, 

y = ^{q*+p*) + x, v y = v. 



In terms of the new variables, 



2 

One checks, 



o(p*-<n-(^+x) 



+ 2 



q* = \{p*-q*)~\v 

P* = \{p*-q*)-\v 

V = -2x 

x = v-\ip*-q*) + \r 



(18) 



(19) 



(20) 



The 3'rd equation reproduces the secondary constraint, and the 4'th equation 
fixes the parameter 77 as expected. 

The above shows, that as a consequence of having the conjugate pairs in 
strong involution, we can simply set the constraints in H equal to zero and 
proceed with the Hamiltonian 



Ho = \(f-q*) 2 - 
Hence we are left with only one HJ equation, 

dS l,dS 



dt + A K dq* 



*\2 



-q*) 



(21) 



(22) 



showing that S will be just a function of q*. The HJ equation associated 
with the primary constraint V = tells us that S is independent of the new 
coordinate x- *P — ff — 0- 

Since the Hamiltonian does not depend explicitely on time, we must have 

dS 2 

—— = —a = constant 
dt 

so that equation (22) has the solution 

S(q*,a,t) = ±2aq* + ]-q* 2 - a 2 t . (23) 

Notice that S depends only on one integration constant, since when regarded 
as the generating function taking us to new momenta p* and V, one of the 
momenta (V) is just the primary constraint V = <fi = 0. There is another 
constant arising from the equation 

da ' 

which implies from (23), 

P = ±2q* - 2at . 

Together with (we choose the plus sign in (23), without loss of generality) 

p* = — = 2a + q* 
oq* 

we now solve the problem in the standard way. Taking account of the con- 
straints Xa = V a = in (18), we find 

1 1 

x — -(3 + at , y = (a + -0) + at 

which evidently is a solution of the Euler-Lagrange equations associated with 
(13). 

Example 2: QM a - model 

Consider now the quantum mechanical version of the "sigma model" as 
given by the Lagrangean 

■J 

L= q - + \{cf- 1) (24) 
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implying the primary constraint p\ = 0. The canonical Hamiltonian on the 
primary surface reads, 

H = ^-X(f-1). (25) 

The complete set of primary and secondary constraints can be written in the 
form 

(f) = p x , ip = X + -, ^1 = ^-1, ip 2 = • (26) 

The Poisson brackets of the constraints show that they can be grouped into 
mutually commuting, canonically conjugate pairs: 

{</?,</>} = 1, M,^ 2 } = 1. (27) 

Let us consider for simplicity the case of two dimensions. One may be 
tempted to make a transformation from (\,p\) to (ip,(f>) and (qi,q2,Pi,P2) 
to (if i, ip 2 , 0, J), where (6>, J) are the canonically conjugate variables, 

= arctan(g 2 /gi) , J = e^qipj , (28) 

with J evidently playing the role of the "third" component of angular mo- 
mentum. The second set of variables commute with the first set, except for 0, 
which does not commute with </?: {(/?, 0} = J/(f. Hence this transformation 
is not canonical. 

In order to remove this problem we eliminate A and p\ using the con- 
straints (p = and <fi = p\ = 0. 7 The resulting Hamiltonian is H = \p 2 q 2 . 
Noting that 

J 2 = ff - (q- pf =f(l + Vl ) - 4^ + l) 2 
we then obtain for H , 



H = ~-zJ 2 



1 

2 

We are thus left to solve the equation 

2 



dS 1 dS\ 
I + 2l =° < 29 > 



7 It is easy to see, that this is a legitimate procedure since A plays the role of a Lagrange 
multiplier in the Lagrangean. More precisely, implementing strongly the constraint ip = 
in the Hamiltonean (25) from the outset, leads to an equivalent set of equations of motion. 
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which has the solution (an irrelevant additive constant has been ommitted) 



a 2 



S(9,a,t) = —-t + a9. (30) 



In order to obtain the solution to the Euler-Lagrange equations we recall 
that we still have the equation §f = which gives the solution in terms of 
two integrations constants: 

d = at + (3 (31) 

One again verifies that this provides the solution of the Euler-Lagrange equa- 
tions. 

Example 3: Multidimensional rotator 

Consider the Lagrangean [8, 11], 

L = \p + Xq- q- 

There is one primary constraint = p\ = and the canonical Hamiltonian 
evaluated on the primary surface reads, 

H Q = l -{p-\<!)\ (32) 

The model effectively describes the motion on a n — 1 dimensional sphere 
without specification of the radius of the sphere. Note that unlike in the case 
of the (T-model, the momentum conjugate to q is no longer the mecanical 
momentum: p = q + Xq. There is only one secondary constraint 

up = X - — , (33) 
q 2 

which has been chosen to be canonically conjugate to 0: {</?, <p} = 1. 

i) We first consider this model in two dimensions. A coordinate change 
(A, pa) (v 9 ) 0) an d (Qi,Pi,Q2,P2) —> ( r ,Pr,@,J) would again lead to the 
second set not commuting with the first, so that A should first be eliminated 
using the constraint tp — 0. Having done this we may replace the canonical 
Hamiltonian (32) by (see footnote 7) 



*-H'-W-£ (34) 



7-2 
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where r = y/q 1 , J is given in (28), and {6, J} = 1. Notice that r just enters 
as a parameter. Replacing J by (§§), we obtain from (34) the HJ equation, 

dS 1 fdS\ 2 n 

w + ^[w) =°- (35) 

which just differs from (29) by r being arbitrary, and not equal to one. 
Correspondingly (30) is replaced by 

a 2 

S(6,a,t) = t + raQ. (36) 

The rest proceeds as in the case of model 2. 

ii) We now reconsider the problem in n-dimensions. Solving again the 
constraint (p — for A, the HJ equation (35) is replaced by, 8 

^ + i f^V _ I = 

dt + 2 ^gj 2 [dqj [dqj q 2 
We make the Ansatz, 

S = f(n-q) (37) 

with n a unit normal vector parametrized by n—1 constants. The HJ equation 
for / then reads, 

Setting x = n ■ q and r 2 = (f, we then find 



/'(*) = ±- 



so that upon integration in x the Hamilton principal function takes the form 

b = artan —j^^=^^= — — t + const (38) 
'r 2 — (n • q) 2 2 



The Hamilton principal function contains n independent constants, which we 
take to be a and ni, n 2 , . . . n n _i, while the normalization of n implies for the 

8 Note that the coordinates q a play the role of <?* in section 1. 
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n'th component, n n = y 1 — Y%=1 n l- Differentiating the Hamilton principal 
function with respect to these constants (new momenta in the corresponding 
generating functional) yields the n time-independent new coordinates: 

<99 (99 
0=7T> a = l,2,...n-l. (39) 

From the first equation and (38) we easily obtain 

n - q = r cos — — °^ = r cos Q(t) . (40) 
Using this, we obtain from the second equation in (39) 

g a = ^sinfi(t) + ^; a = 1, 2, . . . , n - 1 (41) 

Multiplying this equation by n a , summing from a — 1 to a = n — 1 and using 
the normalization of n we obtain from (41), 

n ■ q = — ^ n aPa sin Q(t) H — — —r cos f2(t) 

Q! a=1 n n 

This allows us to solve for q n as 

/ l n -l \ 

q n = n n \r cos V n a/ 9 a sin Q(t) . 

Substituting back in the equation for q a we obtain 

1 ( n-i \ 

q a = - [ (3 a - n a Y^ n a f3 a sin Q(t) + n a r cos Q(t) , a = 1, 2, . . . , n - 1 . 



a 



a=l 



— # 

Introducing the n-dimensional vector /3 = (Pi, fa, ■ ■ ■ , Ai-i,0), these results 
can be written compactly as 

q— — (0—(n- 0)n) sin f2(£) + rn cos Q(t) 

Substituting the above result into our original condition, r 2 = q 2 , leads to 

/3 2 - (ra • /3) 2 



(42) 
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Thus the radius of motion is fixed if the new time independent coordinates 
are specified. More conventionally one would consider the radius of motion 
as a given initial condition; then the above condition eliminates one of the 
P a as an independent constant. Indeed, from (42) we see that the component 

— * 

of P orthogonal to n is fixed to be ar. In this case one therefore considers 
a, r and the remaining n — 2 constants, which determine the component of 
P parallel to n as independent constants. 

Example 4: Landau model in the zero mass limit 

The Lagrangean of a spinless charged particle moving on a two dimen- 
sional plane in a constant background magnetic field B perpendicular to the 
12-plane, and a harmonic oscillator potential, can be written in the form 

m B _ k _q 

Liandau = —q + —q x q - -q . 

In the zero-mass limit the Lagrangian of the model reduces to 

L = ^qx?-±? (43) 

There are only two primary constraints 

A, 1 M ^ 



with non-vanishing Poisson brackets 

{4>i,4>j} = en 

We perform the following change of variables from qi,Pi to 0i, 2 , q,p: 9 



1 1 B 

q = 7> qi + ~b P2 ' V = V\~ , (44) 



3 In the notation of section 1, q,p now play the role of q*,p*, and <j>i — x, 4*2 = V. 
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with the inverse transformation 



Qi = 3 -^02, = ^(VB<Pi - P) , 

Pi = ^(y/Bfa+p), P2 = \(Bq + ^B<p 2 ). (45) 

In terms of the new coordinates the canonical Hamiltonian evaluated on the 
constraint surface reads, 



k ( , 1 
2 V" ' B~ 21 



Ho=~W + ^P 2 ) ( 46 ) 



Hence we have for the HJ equation 



dS k / o 2 



having the solution 



/■? /2q; 

a, t) = -at + B J dq' Vy - g /2 + const. (48) 
From the equation (3 = |^ we then obtain in the usual way the solution 



with uj = k/B and b = f3k/B. 



3 Conclusion 



A nazwe extension of Hamilton- Jacobi (HJ) theory to second class constrained 
systems is condemned to failure right from the outset since it leads to differ- 
ential equations which are in direct conflict. Mathematically this conflict is 
expressed by the violation of the integrability condition (6). 



14 



It is useful to reformulate these integrability conditions in terms of the 
commutator algebra of operators. To this end we begin by introducing for 
each H' a a linear differential operator X„ defined by [6, 9, 10] 

XJ = {H' & , /} , (49) 

where the Poisson bracket is the generalized Poisson bracket defined in (7), 
and / an arbitrary functions of and p^. If H' a does not possess an explicit 
time dependence, we recover the usual definition of the Poisson bracket. Mak- 
ing use of the Jacobi identity for Poisson brackets, we obtain from Eq.(49) 
the following relation between the commutator of the operators and the 
Poisson algebra of the constraints : 

[X*,Xp\f = {{H' B ,H£,f}. (50) 

for arbitrary /. It is now clear that the commutator algebra of the operators 
Xa will not close unless the Poisson algebra of the H' a closes. 

One way in which this problem has been approached is to iteratively 
extend the set of operators H' a until the algebra closes. Equivalently, new 
operators in (50) are iteratively introduced, enhancing the original set X f± to 
a larger set X„ of operators until the closure relation 

[X 5 ,X^]=c^. (51) 

is achieved [12]. The corresponding system of partial differential equations 
will then be integrable. However, there is no general reason why the solution 
for the Hamilton principal function S for this extended first-class algebra 
should provide the solution for the problem in question. 

As an alternative approach to formulate a consistent set of HJ equations, 
one can consider [8] the embedding of the second class system into a first 
class one characterized by a strongly involutive algebra, following the BFT 
[4] construction. This was not the point of view to be taken in this paper. It 
was rather our objective to confront ourselves directly with the problem of 
non-integrability, and to point out a procedure for coping with this problem. 
The procedure consisted in making a canonical transformation to a new set of 
variables in which the second class constraints become part of the variables, 
grouped in canonically conjugate pairs. As we illustrated by a number of 
examples, the new Hamiltonian is then obtained by setting the constraints 



15 



strongly equal to zero, and the HJ equations associated with the constraints 
are no longer in conflict with each other, nor with this new Hamiltonian. 
The solution for the coordinates obtained in this way were shown to coincide 
with the solutions of the corresponding Lagrange equations of motion. 

If the constraints are nonlinear, one is in general not able to perform a 
canonical transformation, turning these constraints into canonical pairs [7], 
and our analysis as presented here is no longer valid. Investigations in this 
directions are presently under way. 
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